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The methods of quantum chemistry and solid state theory to solve the many-body 
problem are reviewed. We start with the definitions of reduced density matrices, 
their properties (contraction sum rules, spectral resolutions, cumulant expansion, 
iV-representability) , and their determining equations (contracted Schrodinger equa- 
tions) and we summarize recent extensions and generalizations of the traditional 
quantum chemical methods, of the density functional theory, and of the quasi-particle 
theory: from finite to extended systems (incremental method), from density to den- 
sity matrix (density matrix functional theory), from weak to strong correlation (dy- 
namical mean field theory), from homogeneous (Kimball-Overhauser approach) to 
inhomogeneous and finite systems. Measures of the correlation strength are dis- 
cussed. The cumulant two-body reduced density matrix proves to be a key quantity. 
Its spectral resolution contains geminals, being possibly the solutions of an approx- 
imate effective two-body equation, and the idea is sketched of how its contraction 
sum rule can be used for a variational treatment. 
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I. INTRODUCTION 



In the electron theory of atoms, molecules, clusters, and solids (N electrons are confined 
to a volume Q) the pairwise Coulomb repulsion between electrons together with the Pauli 
'repulsion' between electrons with parallel spins causes the complex phenomenon of electron 
correlation fl-js^. The problem is, to calculate this electron correlation from first principles. 
Even within the simplifying non-relativistic description and within the Born-Oppenheimer 
approximation, this problem can not be solved exactly. In some cases, correlation beyond 
the Hartree-Fock approximation causes only quantitative changes, but there are other cases, 
where it makes qualitative changes. So, electron correlation let exist the dimers F 2 and Hg 2 , 
or it gives the correct sign for the dipole moment of CO. Electron correlation allows a solid 
to be a metal (in the Hartree-Fock approximation the density of states vanishes at the Fermi 
energy according to ~ 1/ln \e — £f|)- Correlation brings the Curie temperature of ferromag- 
nets down to the experimental values (in the Hartree-Fock approximation they are much 
too high). Electron correlation is important for the ground states of transition- metal oxides 
to be insulating or metallic, ferro- or antiferro-magnetic. Correlation has its peculiarities 
or needs different computational methods for closed-shell systems and open-shell systems, 
for finite and extended systems, for systems near and far from the equilibrium (defined by 
vanishing Hellmann-Feynman forces), for the ground state (GS) and for excited states (for 
a solid, this means its quasi-particle band structure), for different types of solids, for weak 
or moderate correlation (as for sp metals and semiconductors) and strong correlation (as for 
systems with open d-shells or open f-shells, e.g., transition metal oxides and rare earth or 
actinide systems, cf. Appendix). Here we focus on the GS. Its correlation shows up in the 
many-body (MB) wavefunction (WF) of the GS (e.g. in a Jastrow factor in front of a Slater 
determinant) or in the correlation energy (being unfortunately not sensitive to distinguish 
between weak and strong electron correlation) or in quantum-kinematical quantities, which 
follow from the MB WF by integration over most of the variables leaving only a few of them 
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unintegrated. Thus hierarchies of reduced densities (RDs) pi,P2, - '' an d reduced density 
matrices (RDMs) 71, 72, • - • arise with the advantage to allow a distinction between weak and 
strong correlation. Therefore they are used to define measures of the correlation strength. 
One possibility is the correlation induced non-idempotency of the 1-body RDM (1-matrix) 
7i. Another possibility is to derive correlation influenced particle-number fluctuations from 
the pair density (PD) pi in parts of the system; the conclusion is: correlation suppresses 
such fluctuations lj], p. 157, 2^|. 

Important issues of the RDM theory are listed in the following: the contraction proper- 
ties of the RDMs 71, 72, • ' ' an d their cumulant expansions (Sec. 2), parts of the total energy 
as functionals of RDs or RDMs (Sec. 3), the hierarchy of contracted Schrodinger equations 
(CSEs) [what means the equation for 71 contains 72 and 73, the equation for 72 contains 
73 and 74, • ■ •] (Sec. 4), and the size-extensivity [because only size-extensive quantities are 
compatible with the thermodynamic limit (TDL), which means N, Q — > oo,N/Q = const] 
(Sec. 5). Finally a comprehensive summary of recent developments in the calculational treat- 
ment of the many-electron problem is given (Sec. 6). Thereby, "old" methods to calculate the 
electronic structure like the WF based methods [namely, configuration interaction (CI), cou- 
pled cluster (CC), M0ller-Plesset (MP), and quantum Monte Carlo (QMC) for finite systems] 
and the density based density functional theory (DFT) 2^, 25 1 and the Green's function 
based quasi-particle theory (QPT) 2f| for both finite and extended systems are compared 
with some "new" methods like the WF based incremental method 

nearsightedness of electron correlation jjjj is used) or the 1-matrix based density-matrix 
functional theory (DMFT ) (b eing conceptionally simpler and yielding more informations 
than the DFT) t 

Appendix). Within the RDM theory , there seem to be three possibly prospective 

(?) devests: Qj the iV-representability stud.es (Kummer variety, Coleman, algorithm 
linear inequalities) |8[, (ii) the CSE method for finite systems (jj, p. 85, 117, 139 
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29j (thereby the 
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0, , Q| , and also extensions towards strong correlation (Sec.VI(c),(i) 



including critical discussions of the thereby used 'reconstructions' 
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and alterna- 



tive suggestions (Kutzelnigg |4j|), and (iii) the Kimball-Overhauser approach for the PD of 



the homogeneous electron gas (HEG) in terms of 2-body WFs (geminals), cf. 
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and refs. therein. - The paper contains also open questions with the intention to stimulate 
further research. 

TABLE I: List of acronyms 



AGP 


antisymmetrized geminal power 


MB 


many-body 


CC 


coupled cluster 


MP 


M0ller-Plesset 


CI 


configurational interaction 


MR 


multi reference 


CSE 


contracted Schrodinger equation 


PDFT 


pair density functional theory 


DFT 


density functional theory 


QMC 


quantum Monte Carlo 


DMFT 


density matrix functional teory 


QPT 


quasi-particle theory 


ELF 


electron localization function 


RD 


reduced density 


GGA 


generalized gradient approximation 


RDM 


reduced density matrix 


GS 


ground state 


SR 


sum rule 


HEG 


homogeneous electron gas 


TDL 


thermodynamic limit 


KS 


Kohn-Sham 


XC 


exchange-correlation 


LDA 


local density approximation 


WF 


wave function 



II. BASIC QUANTUM KINEMATICS 

The GS of a non-relativistic iV-electron system in Born-Oppenheimer approximation 
is described by a iV-body WF • • • ,N) with 1 = (r^oq) and j dl = j d 3 ri^2 ai or 
equivalent ly by a hierarchy of hermitian RDMs 71, 72, 73, • ■ ■ [first introduced by K. Husimi 



44j ] with the following contraction sum rules (SRs) 

7l(1 |, H /_^_ 72(1 |,, 2 | 2) , 72(1 | O | 20 ./^ 73(lM2 .,3|3),.. (!) 

or C72 = (N — 1)71, C73 = (iV — 2)72, • • • for short, i.e. the contraction operator C makes the 
last variable pair diagonal (i' = 1) and integrates it (J di). This chain 7at(= 
73 — ¥ 72 — ► 7i is called iV-representability, = ty*(V, ■ ■ ■ ,N'). It becomes infinitely long 
in the TDL. Note Tr7i = N, Tr7 2 = N(N — !),-••, Tr7^v = N\. The diagonal elements of 



lit 72, 73, ' ' ' axe the density pi, the pair density (PD) p 2 [containing the Fermi hole for the 
spin-parallel electron pairs and the Coulomb hole for the spin-antiparallel electron pairs], 
the triple density p 3 , ■ • ■ : 

Pl (l)= 7l (l|l), p 2 (l,2)= 72 (l|l,2|2), p 3 (l,2,3)= 7 3(l|l,2|2,3|3),--- (2) 

or pi = Dji, p2 = D7 2 ,p 3 = D73, • ■ • , i.e. the diagonal operator D makes all the variable 



pairs diagonal (l'=l, 2'=2, • ■ •). Note Pi > 0. I 



Daudel loge, refs. on this loge theory are in Ref. 



X is a certain part of the system (e.g. a 



3,Q 



then 



f , 1 \ (N\ f dld2 . . 
N x = j x dl Pl (l), [ 2 ) x = / x ^T^(l,2),--- (3) 

have respectively the meaning of the average number of particles in X, the average number 
of particle pairs in X etc., approaching N, ( 2 ), • • • for X — > 00. The term 'average number 
of particles in X' means: the number of particles in X fluctuates. In lowest order, such 



fluctuations are described by the mean squared deviation ANx- They are determined 



the 2-body function w 2 (l, 2) [defined by P 2(l, 2) = pi(l)pi(2) — ^(l, 2)] according to [23 

^^ = l-^j x dljd2w 2 (l,2), J dl J d2 w 2 (l,2) = N. (4) 

w 2 (l,2) is size-extensively normalized. Particle-number fluctuations can be described in 
more detail by Px(M), the probability of finding M{= 0, 1, • • • , N) particles in X with 

N N N 

J2Px{M) = l, )T Px(M)M v v- ^Px(M)(M-iV x ) 2 = (AiV x ) 2 . (5) 

M=0 M=0 M=0 

The example j^] 

Px(l)= [ dlpx(l)-i I dl I d2p 2 (l,2) + i f dl f d2 f d3p 3 (l,2,3) (6) 

shows, how e.g. Px{l) follows from the reduced densities. For other values of M cf. j^. 
An approximate expression needs only p\ and p 2 , namely Px{M) ~ e -/3(M-«) 2 w [Hy ^j s _ 
tribution parameters a and f3 determined by Nx and ANx- [This expression maximizes 
the entropy — ^2 M Px(M) In Px(M) under the constraints of given Nx and ANx-] The 
fluctuation analysis of finite and extended systems based on this approximation shows that 



strong) correlation (strongly) suppresses fluctuations' l|, |23|; for attractive interaction (of 
the Calogero-Sutherland model = one- dimensional Fermi gas with 1/xfj interaction), they 



are enhanced |47[ . For the peculiarities of strong electron correlations cf . Appendix. - The 
'on-top' curvature of the Fermi hole is used as a local measure of the correlation strength 
It is also derived an 'electron localization function (ELF)' 48^ an alternative is introduced 
in Ref. 3- 

Because the RDMs are hermitian, they can be diagonalized with 



7i 



N 



(7) 



as the simplest example. The ip K (l) are called natural orbitals (more precisely spin-orbitals) 
with natural occupancies u K . For the HEG, the ip K (l) are plane waves and the v K yield the 

occupancies are only 

exceptionally idempotent, in general they are non-idempotent: v\ < v K and ^ K v 2 K < N. 



Because P.-O. Lowdin has asked for the meaning of Y2 K v \ = 



511 ]. the quantity 



C-2 



< c 2 < 1 



(8) 



will be referred to as Lowdin parameter (with tr = ^pTr). For non- interacting single- reference 
cases (only one Slater determinant) and therefore idempotent occupancies v K = or 1, 
this parameter vanishes. Interaction then causes a correlation tail in the CI expansion 
with a lot of additional Slater determinants beyond the leading Slater determinant. This 
makes the v K non-idempotent, thus c 2 > 0. So, c 2 is another measure of the correlation 



strength 
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53| . An increasing weight of the CI correlation tail increases also c 2 (with 
perhaps weak, moderate, and strong correlation for c 2 -C 1/2, c 2 ~ 1/2, and c 2 ^> 1/2, 
respectively). Thereby the CI expansion is thought to be unique (basis set independent) by 
representing it through natural spin orbitals. Often one has to start with a multi-reference 
(MR) state, a sum of several Slater determinants to correctly describe, e.g., a dissociation 
limit or nearly degeneracy. Then the Lowdin parameter is already non- zero for 'no CI 
correlation tail' (cmr > 0. For this descript ion t he term 'non-dynamical or static correlation' 
is used) following Clementi and Corongiu lC^ . Full CI with its additional correlation tail 



('dynamical correlation') makes then C2 > cmr- - For a solid, v K contains the occupancy 
band structure 5^ |. - The 2-body analog of Eq.(|7]l is the spectral resolution of the 2-matrix 
72, containing antisymmetric 2-body WFs (natural geminals). Are they the solutions of a 
2-body Schrodinger equation with an appropriately screened Coulomb repulsion (similar as 
this is the case with the Kimball-Overhauser approach for the HEG 411])? 



An important aspect of the quantum kinematics is the cumulant expansion of RDMs and 



RDs |55j. Its most compact representation uses generating functionals, cf. e.g. P. Ziesche 
in Q], p. 33. For the simplest case of the 2-matrix this means 

72 (l|l',2|2')= 7 i(l|l / )7i(2|2')-7i(l|2')7i(2|l / )-X2(l|l / ,2|2') (9) 

or 7 2 = A7171 — Xi f° r short (A = antisymmetrizer). The factorized (or disconnected) terms 
are the generalized (because the natural occupancies are non-idempotent) Hartree-Fock part 
of 72 with Tr 7 2^ F = iV(iV — 1) + Nc2- The non-factorizable remainder \2 is called cumulant 
2-matrix. In perturbation theory, X2 is given by (size-extensive) linked diagrams; therefore 
it is also called connected part of 72. [An equation similar to Eq.(JH} defines the cumulant 
part of the 2-body Green's function.] The above mentioned contraction SR C72 = (iV — l) 7l 
(which is not size-extensive) is transformed to C%2 = 7i — Cjf or 

J d2 X2 (l|l',2|2) = 71(1110- / rf2 7 i(l|2) 7 i(2|lO = ^^(l)^(l-^)C(lO- (10) 

The invariance of the rhs under the exchange v K (1 — u K ) is called particle-hole symmetry 
56^ . Eq. fjlOj) contains the normalization SR 

tr%2 = ^JZ^(1 - ^k) = c 2 ; s = --^^[z/„lni/ lt + (1 - ln(l - v K )\. (11) 

K K 

Thus the Lowdin parameter C2 measures the non-idempotency of the natural occupancies u K 
and simultaneously the normalization of the cumulant 2-matrix \2 per particle, thereby not 
affecting the norma 
correlation entropy 



ization of the total 2-matrix 72 . s defines the particle-hole symmetric 
5(1 152I ] . For the PD P2 the cumulant expansion Eq. (jHJ) means 



p 2 (l,2) =p 1 (l)p 1 (2)- 7l (l|2) 7l (2|l)- U2 (l,2), w 2 (l,2) = X 2(l|l,2|2). (12) 



m 2 may be addressed as cumulant PD. Comparison with Eq.(JU) shows u> 2 (l,2) = 
7i(l|2)7i(2|l) +1*2(1)2), being thus the exchange including cumulant PD. 

If Xi is available from perturbation theory or from some effective 2-body theory, then the 
SR Eq. fjlOjl allows one to determine the 1-matrix 71 by means of diagonalizing Eq. fjlOjl . The 



resulting quadratic equation A K = z/ K (l — u K ) is solved by = \ ± y \ — X K . Associating 
certain k's to + and all the others to — , yields occupancies v K > | and v K < | respectively, 
such that XIk^ = N. [For the low-density HEG(r s ^> 1), i.e. for the Wigner crystal only 
v~ has a physical meaning and it is X K -C 1, thus v~ ps X k .] Let us denote the described 
procedure by 71 = C%2- With this 71 the exact MB kinetic energy T and the Lowdin 
parameter c 2 can be calculated and also the PD p 2 , which is needed for the interaction 
energy V\ n t and for the fluctuation analysis (both in direct and momentum space). So, the 
cumulant 2-matrix \i is f ne most interesting quantum kinematic quantity. Although it is a 
key quantity, unfortunately not much is known about it. Questions: 

• What is known about the spectral resolution of the cumulant 2-matrix % 2 ? 

Xa (l|l', 2|2') = I ^jc(1. 2W«1', 2'), TV E ^ = c 2- ( 13 ) 

K K 

Which values are possible for the cumulant occupancies /x^? 

• 7 2 is A^-representable Jg|. What are the consequences for % 2 ? 



What are the peculiarities of X2 for metals (for the high-density HEG cf. |57|), quasi- 
one-dimensional charge- density- wave conductors, semiconductors, (band, Peierls, 
charge-transfer, Mott-Hubbard) insulators, superconductors (off-diagonal-long-range 
order), ferro/antiferro- magnets, ferroelectrica? For cases of strong electron correlation 
cf. Appendix. 



What jire the peculiarities of \2 for cases, where one of the occupancies of 7 2 is very 
large 
cf. 



£ [.There exists a rich literature on antisymmetrized geminal power (AGP) WFs, 
I. M. Rosina has studied the 2-matrix 72 of the AGP WF ml. The AGP 

n 

function is a flexible ansatz for fermion systems with arbitrary [8]. Extreme AGP 



1, p. 35, Q 



WFs have (for large N) very large 2-matrix occupancies equaling [N/2] 
The GS of Be has a large 2-matrix occupancy (> 1) (V.H. Smith, priv. commun.). 
What does all this mean for the cumulant occupancies p^-? 



Electron correlation is discussed in terms of its nearsightedness 
show up by % 2 ? 



[30]. 



How does this 



A drawback of X2 is that the positivity condition p 2 > becomes more complicated. 

The normalization and contraction of X3 yield Tr%3 = 2^ K z/ K (l — z/ K )(l — 2u K ) = NC3 
and Cx 3 (l|l',2|2',3|3') = 2x 2 (l|l', 2|2') - [A'Cx 2 (l|l', 2|3)xi(3'|2 / ) + h.c], where the con- 
traction operator C makes 3' = 3 and J d3 as defined above, A' is the antisymetrizer with 
respect to the primed variables and h.c. means hermitian conjugate. Thereby integrals 
/ d3^(l, 3)^*(3) appear, which describe the overlap of a natural orbital with a cumulant 
geminal. 

III. BASIC FUNCTION ALS 

Explicitly defined functionals for the components of the GS energy E are 

T[yi] = J ^(rxhiail')!^, *(r x ): 



1 d 2 



2 • 



2dr\ 

dl /^(lHxtOx), V int [p 2 ] = / -y- p 2 {l,2)v int {r 12 ). (14) 

These are linear functionals. The generalized Hartree and Fock components of the interaction 
energy V- int are the bilinear functionals 

Vu[pi] = J Pi(l) P i(2)v int (r 12 ), 7 P [ 7 i] = / 7i(l|2)7i(2|lKt(ri 2 ). (15) 

Furthermore, the cumulant partitioning needs the linear functionals 

VfcM = - / ^r w 2(l,2Kt(r 12 ), Vb[u 2 ] = - / ^ « 2 (l,2Kt(r 12 ). (16) 



Then the interaction energy is given by 



fd2p 2 , 



MntN = Vk[pi] + V FC [w 2 }, pi = — _ - , w 2 = p l - p 2 (17) 



or 

Mnt[7i, u 2 ] = V H [pi] + V F [jx] + V c [u 2 ], Pl = D 7l = L^l. (i 8 ) 
Consequently, the total GS energy is 

£[ 72 ] = T M + KxtN + Vkl>i] + ^fcM (19) 

or 

% 2 ] = T[ 7l ] + KxtN + ^hN + vm + V c [u 2 ]. (20) 

The N 2 terms of Kxt and Vh cancel each other. In the sum Vh + Vp, the self- interaction 
terms cancel each other. In the sum Vp + Vq there is also a certain cancellation as it may be 
seen from the normalization Tr(7 2 + W2) = N, which does not contain the Lowdin parameter 
C2. To find the GS, one has to minimize these functionals E[ r j 2 ] or £"[x2], provided that the 
iV-representability is taken into account: 72 follows from 73, 73 from 74 etc. by successive 
contractions, cf. Eq.JT]). 

n 

Other functionals implicitly defined by means of Levy's constrained search [2QJ| are (i) for 
the well tried and widely used DFT 

(T + V int )[ Pl }= wm(*\f + VW|*), Pi = D 7l , 7l = -^ 2 -,---, (21) 

from which follows (T + V FC ) [pi] = (T + V int ) [pi] - V n [px], (ii) for the DMFT (under dis- 
cussion) 

^[Ti]=ndn(*|y int |*), 7l = -^-,---, (22) 

*^7l iV — 1 

from which follows Vc[7i] = Vint [71] ~ Vk[pi] — ^f[7i], and (iii) for a possible pair-density 
functional theory (PDFT) 

f[p 2 ] = min(*|T|*), p 2 = D 7 2, l2 = -^-,---. (23) 

The 1-matrix functional Vint [71] has the simple scaling property Vi nt [7i] = AVj nt [7i], where 
7x( r i| r 'i) = A 3 7 1 (Ar 1 |Ar' 1 ). The PD functional T[p 2 ] has the simple scaling property T\p%\ = 
A 2 T[p 2 ], where P2( r i5 r 2) = A 6 p 2 (Ari, Ar 2 ). The reason of these simple scaling properties 
are the homogeneities of the corresponding operators V{ nt and T with degrees —1 and —2, 



respectively. These different degrees cause the more complicated scaling of the 



functional (T + Vi nt )[pi], where the constrained search involves both T and V int 61]. To 



density 



one may search for approximate 
and 7i[p2], respectively. From the 



find approximate functionals for Vint[7i] and for T\ 
functional X2I71] (as e.g. proposed by K. Yasuda [5_ 
above introduced functionals follow the density functional _Edft[pi] = (^+^4nt)[pi] + V r C xt[pi], 
the density-matrix (1-matrix) functional £t>mft[7i] = ^[71] + Kxt[pi] + ^int[7i], and the PD 
functional £pdft[P2] = T[p 2 ) + KxtL°i] + ^inttpa], to be minimized, cf. Sec.VI(c),(d),(e). 



IV. BASIC EQUATIONS 



In order to find the GS \T/ ( 1 , • • ■ , N), the starting point is the MB Schrodinger equation 



3 



by multiplying it with 



= E^. ^From it follows the hierarchy of CSEs 
\E'*(1 / , • • • ,N') and performing successively the N — 1 contraction, the N — 2 contraction 
etc.: 

d2 ■ ■ ■ (H - E) **'*| 2 , =2i3 , =3! ... = 0, Jd3---(H-E) **'*| 3 , = 3 !4 , =4i ... = 0, • • • . (24) 

On the lhss the RDMs 71, 72, 73, • ■ • appear. Thereby the CSE for 71 contains 72 and 73 (1st 
formula), therefore it is called 1,3-CSE; the CSE for 72 contains 73 and 74 (2nd formula), 
called 2,4-CSE, • • • . Expressing these RDMs by their cumulants X2, X3, • • • an equivalent 
hierarchy appears (no longer containing E); for a uniform notation one should define the 
1-body 'cumulant' as Xi = 7i- 

An equivalent writing of the Schrodinger equation is E = min(^lHl^), (^l 1 ! 7 ) = N\. An 



equivalent writing of the hierarchy (|2"4"j) is 

E = mmE[ l2 ], 7a = -?2L, . . . , Tr 7i = N(N - 1) • • • (N + 1 - i). (25) 
72 A — z 

The size-extensive version of Eq. (j23j) is e = mine[x 2 ]) X2 = Cx3, • • • , trji = Ci with 

X2 

e[ X2 ] = E\x2\/N, cf. Eq. pnjl . It survives the TDL (N, Q -> 00, N/tt = const). 

Using the relation ji = CX2 described after Eq. flHj) , the total energy becomes a functional 
of X2, which in its spectral resolution provides geminals tpx and occupancies fix as variational 
parameters. Does this procedure give a link to the Kimball- Overhauser approach? 



V. SIZE-EXTENSITIVITY AND THERMODYNAMIC LIMIT 



Solid state theory considers extended systems, which result from finite systems (N elec- 
trons confined in a volume Q) through the TDL. In the terms of thermodynamics one has 
to distinguish intensive quantities, e.g. the density pi(l), and extensive quantities, which 
are proportional to N, like the normalization of p\ or the kinetic energy T. Amongst the 
quantities listed above are also such ones, which contain terms proportional to N 2 . Exam- 
ples are V ext and V^. Fortunately their sum is size-extensive, what makes the total energy 
E also size-extensive. Thus the ratio e = E/N survives the TDL, it defines the bulk energy. 
Another example is the normalization of the PD p 2 , where also N 2 terms are present. This 
deficiency is 'removed' by considering the cumulants of p-i. Both u>2(l,2) and 7x2(1, 2) are 
size-extensively normalized and the cumulant 2-matrix \2 is size-extensively contractable. 

What happens with the equation H ^ = E^> in the TDL? One may divide it by N and 
use the bulk energy e = E/N on the rhs, but H/N has no meaning in the TDL. Fortunately 
the hierarchy ()24j) [being finite for finite N, Q] only becomes infinitely long in the TDL. 
This hierarchy has to be solved at least approximately. Even if this has been done yielding 
the GS properties X2 ( an d from it also 7i,pi,p 2 ) and e, the question arises, how to obtain 
such important solid state properties as the quasi-particle band structure with its band gap 
and particle-hole bound states (in the case of a semiconductor) and its Fermi surface and 
plasmons (in the case of metals) and its effective masses and life-times. This task goes 
beyond RDM theory. It needs the Green's function based QPT, cf. Sec.VI(i). 

The TDL let emerge solid state properties having no counterparts in small clusters. Such 
properties are long-range ordering phenomena like ferromagnetism or other strong-electron- 
correlation phenomena like metal-insulator transitions, Wigner-charge ordering, coexistence 
of ferromagnetism and superconductivity, heavy fermions, quantum criticality, cf. Appendix. 

VI. CLASSIFICATION OF CALCULATIONAL METHODS 

Amongst the calculational methods one may distinguish methods for finite systems 
(atoms, molecules, clusters) and for extended systems (crystalline solids) including jellium 



models. Another classification is: methods based on the many-body WF and non-WF based 
methods, cf. Table ITT1 



TABLE II: Methods to calculate the electronic structure. The symbols (!), (?), (??) mean 'well 
tried', 'option', 'idea', respectively. 





finite 


extended 




WF 


CI,CC,MP(!) 


? 


(a) 


based 


QMC(!) 


? 


(b) 








(c) 




DMFT(?) 


DMFT(?) 


(d) 


non-WF 


PDFT(??) 


PDFT(??) 


(e) 


based 


Coleman(?) 


_> ? 


(f) 




CSE(?) 


? 


(g) 




? <_ 


Kimball-Overhauser 


(h) 




QPT(!) 


QPT(!) 


(i) 



(a) There are recent attempts to make the well-tried quantum-chemical methods, e.g. 
CC, applicable also to crystalline solids. This subject was thoroughly discussed at 
2(3 ]. especially by P. Fulde, B. Paulus, and U. Birkenheuer (incremental method). 
With this method (based on localized Hartree-Fock orbitals and a 'constrained' CI 



271 ] . One extension of 



scheme), GS properties of insulating solids were calculated 
the incremental scheme concerns metals and strongly correlated systems. Thereby 
the problems appear of how to construct local orbitals in metals, how the incremental 
scheme works with these orbitals and how to handle the high degeneracy near the Fermi 
edges (with a MR method) j&j]. This allows the method to be applied to the strongly 
correlated dissociation limit and describe a metal-insulator transition 2]. Another 
extension of the incremental scheme aims to compute correlated wave functions and 
energies of both valence and conduction bands 




29( . Thus the band gap for periodic 



insulating systems including the effect of electron correlation can be determined. First 
results are available for diamond and polyacetylene y. 

R.J. Bartlett says "CC theory offers the natural vehicle for transferring lst-principle 
electronic structure informations into material modeling", in f3b|, p. 219. Because 

n 

electron correlation is 'short-sighted' (W. Kohn 1301) a localized correlation treatment 
is developed using natural bond orbitals 66J. 

(b) For the QMC method in general cf . Q| . In order to reliably predict materials proper- 
ties, a QMC method using Hubbard-Stratonovich auxiliary fields has been developed 



671 ]. This MB method is aimed at treating electron correlations in real materials. A 
first application concerns bulk Si. Slater determinants from DFT calculations were 
used as the trial WFs. The cohesive energy is comparable to (or better than) the best 
existing theoretical results. 

In contrast to the WF based methods (a) and (b), the methods (c)-(i) are motivated by 
asking "Is the WF not a monster containing too much information, why one should struggle 
with it? Let us better concentrate on the lowest order RDMs, namely 71,72 and their 
positivity properties and more general their iV-represent ability (Coulson's challenge)", ||. 

(c) DFT j^J uses the density p\ [which has its charme and power 2^ ] as basic quan- 
tity. In its Kohn-Sham (KS) version, the non-idempotent 1-matrix 71 is the sum of 
an idempotent part 7f s , corresponding to a non-interacting system with the kinetic 
energy T K s[pi], and a remainder, such that (i) D7f s = ^° c ° |?/> K (1)| 2 (these KS or- 
bitals ip K are not natural orbitals!) gives the correct GS density p% of the interacting 
system (although •yf s is not the correct 1-matrix), (ii) the equation (T + V|nt)[pi] = 
Tks[pi] + Exc[pi] defines the exchange-correlation (XC) functional £xc[pi] to be ap- 
proximated (LDA, GGA, • ■ ■ ), and (iii) the minimization of £tdft[pi] yields an effective 
1-body Schrodinger equation with an XC potential uxc(l) = ^xc[pi]/^Pi(l)- Finally, 
its solution gives (together with the aufbau principle) the density pi(l) and the total 
energy E. Besides Lagrange multipliers e K appear, which surprisingly approximate the 



quasi-particle band structure of sp metals as (from a conceptional point of view) ille- 
gitimate children of the KS-DFT (W. Kohn). But for semiconductors their KS-DFT 
band gaps are too small and for heavy-fermion systems the experimental effective 
masses are much larger than the KS-DFT ones. DFT goes beyond the correlation 
neglecting Hartree-Fock approximation. It is well-tried and widely used for both finite 
and extended (metallic and semiconducting/insulating) systems, cf. e.g. 0. The 
interconnections between DFT and CC theory are studied (R.J. Bartlett at 18j]). Ac- 
tivities are going on trying to remove deficiencies and failures of this method (e.g. the 
van-der-Waals interaction Q], the problems listed under (d), the asymptotic behavior 
of the XC potential j^), and to generalize it, e.g. DFT for MR cases (A. Savin at 
1^1), time-dependent DFT jzi| and its relation to QPT jljj, DFT for excited states 
72j ]. DFT for transport properties 7^. In particular the originally used LDA [based 
on the HEG and being the Thomas-Fermi approximation on a higher level] is modi- 
fied (i) by climbing the 'Jacob's ladder' via GGA and Meta-GGA up to Hyper-GGA, 
which are non-empirical approximations for i?xc[pi] designed for molecules and solids 
jz^j], (ii) by studying orbital functionals (optimized effective potential method, the DF 
method in 3rd generation), cf. e.g. [3], and (iii) by taking into account peculiarities 
of strong correlation (Appendix), which suppresses charge fluctuations, thus localizes 
electrons and therefore favors insulating behavior. Strong electron correlation is often 
treated with the help of the Hubbard model (containing the on-site Coulomb repul- 
sion U), the t-J model, the Anderson model (including attempts to determine model 
parameters from 1st principles 76]) or the orbital model [77( and the density-matrix 



renormalization group |7Js|. One way for (iii) is the combination with the Hubbard 
model (LDA+U); its possible link to 1-body Green's function theory is discussed in 

Ref.Q, 

where a density functional application to strongly correlated electron sys- 
tems is presented. Another way is the combination with the dynamical mean field 
theory of strongly correlated fermion systems JsiJ • This allows realistic calculations of 
transition metal oxides and f-electron materials, e.g. the Mott-Hubbard transition of 
V 2 3 , the volume collapse transition of Ce, the magnetic properties of LiV 2 4 , and 



the photoemission spectra of SrV0 3 and CaV0 3 



(d) DMFT uses the 1-matrix 71 as basic quantity |3Jj. Unlike DFT, it does not need a 
(non-interacting) reference system, but it needs an approximation for Vi n t[7i] or Vc[ti], 
from which follows the functional -Edmft[7i]- Its minimization yields 71 and E (but 
not p2 and in the case of solids not the band structure), cf. e.g. |£2f. Each of the so far 
proposed approximations have their successes and shortcomings (e.g. violation of the 



positivity condition or o 
functionals are in Ref. 



the contraction SR). Conditions for approximate 1-matrix 
32| . From a conceptional point of view, DMFT is simpler than 
DFT. The hope is that the new method is like DFT cheaper than CI, CC, MP and 
applies to both finite and extended (metallic and insulating) systems and that it can 



handle problems where DFT fails 



8;: 



negative atomic ions are not bound, ionization potentials are too small, 

ation energies of He(Z) and Be(Z) have not the correct large-Z behavior, 



corre 
cf. 



• errors in the bond lengths and bond angles of some molecules, e.g. 
CH 4 , H 2 0, BH, O3, Li 2 , Be 2 (through this series, static correlation is progressively 
increasing) , 

• the reaction H 2 + H — > H + H 2 is another MR case, it has no transition state in 
LDA and one that is too low with better functionals beyond LDA, 

• wrong dissociation limits, 

• cohesive energies of bulk metals are overestimated in LDA and underestimated 
in GGA, 

• DFT poorly performs for polyacetylene fragments and oligoporphyrins, incor- 
rectly predicting a triplet GS for polyacetylene Ssj ] - 

There are links between DMFT and the RDM theory (reconstruction of the 2-matrix) 
j^j] and between DMFT and the theory of antisymmetrized products of strongly or- 



thogonal geminals 



For the geminal theory in general cf. 



(e) The possible PDFT uses the PD p 2 as basic quantity. Again unlike DFT, it does not 
need a (non-interacting) reference system, but it needs an approximation for T[p2], 
from which follows the functional -Epdft [P2] ■ Its minimization together with an ansatz 
for p2 in terms of geminals should yield p 2 and E for both finite and extended systems 



87| . - Another method, 



(but not 7x and in the case of solids not the band structure) 
which uses also two-particle states (what involves the extension of LDA from densities 
to PDs), is in Ref.[3b], p. 325 and refs. therein. 



(f) For the RDM theory in general cf. p, |y, LD, |8|. For the investigation of the N- 
representability in terms of Kummer variety, Coleman's algorithm, and linear inequal- 
ities cf. Ref.j8j. Whether this way is really cheaper than a full CI calculation and 
whether it applies to extended systems has to be checked. 

(g) The generation of higher-order RDMs from 72 via 'reconstruction' has been developed 
by Nakatsuji, Valdemoro, and Mazziotti Q|, p. 85, 117, and 139, respectively, and 
refs. therein, cf. also j^, Q, Hi- The philosophy behind this approach is: Let 
us use an exact functional and violate slightly iV-representability rather than using 
approximate DFT functionals. One may criticize this attitude with the argument: 
Already a slight violation of the iV-representability opens regions of the variational 
space, where the calculated energy is below the GS energy [38]. There seem to be 
encouraging results: using the positivity condition on the 2-matrix (as a necessary 
condition of the iV-representability), energies of closed- and open-shell atoms and 
molecules were calculated variationally 3(| and energies of atoms and molecules (both 
near and far from equilibrium geometries) were calculated solving the 1,3-CSE 
But critical discussions concern 'how the results change, when going from minimal to 
extended basis sets?' and 'are the 3- and/or 4-cumulants really negligible? 

Whether the CSE method (provided that it survives these criticisms) applies also to 
extended systems has to be checked. 



(h) In the Kimball-Overhauser approach 



4l| the singlet PD g+(r) and triplet PD g~(r) 



of the HEG (here normalized as g±(oo) = 1) are parametrized in terms of geminals 



Ri(r, k), which are intuitively solutions of 



1 d 2 



-J + 



1(1 + 1) 



ip it 
such that g(r) = \ [g+(r) + 3g_(r)] with 



+ -+< CT (r 
r 



Ri{r,k) = 0, 



(26) 



g ± (r)=2j2^J2^ R "^ k ^ = 4l>(*iMk). (27) 



L k 



K 



where k\ = ||K + k|,fc 2 = |K — k| , L — (I, mi), + is for even / and — for odd 
I. v^. r (r) describes the effective screening of the Coulomb repulsion by the Fermi- 
Coulomb hole around each electron; on the Hartree level it follows from Au* r (r) = 
4?rp[l — q( r)]. The Ri are scattering states, characterized by scattering phase shifts. 
In Ref . |43| the normalization and contraction SRs of the 2- matrix is related to these 
phase shifts. According to these new SRs it should be possible to calculate not only 
the PD from f scr (r), but also the momentum distribution n(k). For further details 



cf. 



41 



42 



and refs. therein. - Can Eq. (j26|) be derived from MB theory (CSE, 
contraction SRs, Bethe-Salpeter equation or other methods)? - To what extent is 
the Kimball-Overhauser approach related to the PDFT (e) and to the 2-body cluster 
expansion of Ref. 8J|? - Whether this method correspondingly modified applies also 
to inhomogeneous and also to finite systems has to be studied. - A general remark of 
E.R. Davidson is "Scientists have not yet learned to think in terms of ( 2 ) geminals 
rather than in A" orbitals" M, p. 97. 



(i) The Green's function based QPT (like the KS-DFT) rests upon an effective 1-body 
Schrodinger equation with the local, real and energy independent XC potential vxc(l) 
of the KS equation replaced by the non-local, complex and energy dependent XC 
part £xc(l!l',£) °f the self-energy 2fil |. There holds an integral relation between t>xc 
and Sxc |89(, Eq.(ll). QPT needs an approximation for the S X c (e.g. the GW 
approximation and its vertex correction) and yields the GS properties 71 and E and 
besides the quasi-particle band structure e t . Total euergies were calculated (EL, 
85) with the Galitskii-Migdal formula 2l| and with the Luttinger-Ward formula 9J,|. 
The real part of Sxc compensates (in the case of semiconductors) for the otherwise 



(namely within DFT) underestimated band gaps and yields mass renormalizations. 
The imaginary part of £xc accounts for life-time effects. The Coulomb interaction 
between electrons in the conduction band and holes in the valence band are included by 
solving the Bethe-Salpeter equation for the electron-hole two-particle Green's function 
9^ |. This provides the excitonic binding energies (9^. The combination with the afore 
already mentioned dynamical mean field theory allows the first-principle approach to 
strongly correlated systems (Sec. VI(c), Appendix) 9^, 95^ . 



VII. SUMMARY 



Summarizing and focusing on solids, there are several methods to calculate such GS 
properties as E, pi, 71, but from a conceptional point of view only QPT Sec.VI(i) and the 
incremental method VI(a) yield also the band structure including Fermi surfaces (of met- 
als) and band gaps (of semiconductors) and effective masses. So far a weak point is the 
PD P2, cf. Sec.VI(e)-(h). This weakness can be removed, if the cumulant 2-matrix \2 is 
calculated. This is a key quantity because from it follows: the 1-matrix 71 (which contains 
the density pi) and using this result, the 2-matrix (which contains the PD). This PD is 
needed to calculate the interaction energy and particle number fluctuations in parts of the 
system (the magnitude of which measures the correlation strength). The spectral resolution 
of X2 leads to geminals, which may be the solutions of a 2-body Schrodinger equation with 
a screened Coulomb interaction if not exact, but possibly to a good approximation. What 
peculiarities has this screening potential for different types of solids like sp metals, transi- 
tion metals, semiconductors/insulators, superconductors, ferro/antiferro- magnets, strongly 
correlated systems? Further work is called for to answer the questions posed in this paper. 
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APPENDIX: CORRELATION STRENGTH AND STRONG CORRELATION 

The strength of correlation is not a matter of the correlation energy. More relevant is the 
quantum-kinematics expressed in terms RDs and RDMs. So the Lowdin parameter c 2 is one 
correlation index, cf. Eq.fjHJ). Electron- number fluctuations ANx in a part X of the system 
are another index, decreasing with increasing correlation, i.e. correlation suppresses such 
particle-number fluctuations HQ , it thus favors electron localization and non-metallic 
behavior. The Calogero-Sutherland model shows that ANx is more sensitive, because it 
distinguishes between attractive and repulsive interaction, whereas c 2 shows practically no 
change when changing the sign of the interaction (the non-idempotency of the momentum 
distribution is qualitatively the same for these two cases) (47 1. 

Cases of strong correlation have to be described by MR states, cf. the discussion after 
Eq.(|HJ). This strong-correlation case is associated with a strong suppression of fluctuations 
or equivalently with strong localization, cf. discussion after Eq.fjHJ). Examples are: (i) 
He(Z) for Z^Z CT = 0.91, where one electron is localized at the nucleus while the other 



one becomes unbound and escapes 



3, (") 



(ii) the H 2 (i?) molecule in the dissociation limit 
R — > 00, where at each nucleus one electron is localized such that the total spin state is 
a singlet (iii) the Wigner crystallization of the low-density electron gas. For the case 



of 'strictly correlated' systems cf. 



96j]. Ferro- and antiferro-magnetism and corresponding 
quantum phase transitions are 'old' strong-correlation phenomena. 'More new' examples 
are: 

• Mott insulators, where the hopping between the lattice sites is suppressed by a large 



Hubbard U: CuO, NiO, CoO |97|, FeO, MnO, CaCu0 2 , V 2 3 , Ca 2 Ru0 4 , La 2 Cu0 4 . 
Examples for correlation-driven metal-insulator transitions are FeSi (under doping), 
LaCo0 3 (under temperature), V 2 3 and Ca 2 Ru0 4 (under pressure). 

Charge ordering in Yb 4 As 3 or NaV 2 5 , charge-density-wave transition in CuV 2 S 4 , 



spin-density waves in Ce(Rui_ :c Rh :r )2Si 2 , charge-spin separation in SiCu0 2 , Yb 4 As 3 
(the thereby appearing low-lying excitations are not quasiparticles, but collective spin 
and charge density fluctuations). 

Quasi-one-dimensional systems: charge-spin separation in polyacetylen (solitons) or 



in Bechgaard salts (heat transport) 
liquid behavior in organic conductors. 



spin-Peierls transition in CuGe0 3 , non- Fermi 



High T c superconductors and other unconventional superconductors [Sr 2 Ru0 4 , 
ScCu(B0 3 ) 2 ]. UGe2, URhGe are superconducting ferromagnets j^]. 

Below characteristic temperatures heavy-fermion systems (mostly compounds of light 
lanthanoides and actinides with their 4f and 5f electrons) show anomalously large elec- 
tronic heat capacities, which can be described as heavily renormalized quasiparticles. 
They have huge effective masses being several hundred times the free-electron mass. 
Examples are (P. Fulde, in [3b], p. Ill): CeAl 3 , CeCu 2 Si 2 , and CeCoIn 5 (Kondo 
effect), Nd 2 - x Ce, E Cu04 (Zeeman effect), CeCue- x An x ; 4f holes (respectively spin exci- 
tations) are decisive for Yb 4 As 3 and 3d electrons for YMn 2 , LiV 2 4 . The compounds 
UPd 2 Al3 and UPt3 are characterized by a combined coherence of localized and itin- 
erant electrons [the 5f electrons are partially localized as in the sandwich molecule 
U(CsH 8 ) 2 ], they are described by a dual model of localized and hybridized 5f elec- 
trons. The compounds CeCu 2 Si 2 , UBe i3 , UPd 2 Al 3 , UPt 3 , Celn 3 are superconductors, 



ZrZn 2 is a superconducting ferromagnet 



mnj . The compounds CeMIn 5 (M=Co, Rh, 



Ir) are antiferromagnets. The compounds YbRh 2 Si 2 and YbRh 2 (Si .95Ge .o5)2 have 
small antiferromagnetic ordering temperatures at 70 m K and 20 mK, respectively, 



and show other peculiar low-temperature behavior 



101]. 



• Quantum criticality near T = (where quantum fluctuations predominate thermal 
fluctuations) in CeCu6-a;Au x , YMn 2 , CePd 2 Si 2 , CoCu^b-zGea,. 

Other aspects are: mixed valence compounds (SmS, Yb4As 3 ), geometrical frustration 
(LiV 2 4 ), frustrated spin systems, non- Fermi liquid behavior (MnSi, CeCu 2 Si 2 , UPd 2 Al 3 , 



YbRh 2 Si2, SrRu03 is a non- Fermi liquid ferromagnet), fractionally charged excitations in the 
fractional quantum Hall effect with its composite fermions (electrons+flux lines), 'fractional 
charges' (topological excitations of special low-temperature phases, possibly in LiV 2 4 ). For 
reviews on strong electron correlation cf. Ref . \v\ . Theoretical treatments (ab-initio or with 
models) are mentioned in Sec.VI(c),(i). 

For solids, the bandwidth W is compared with the Hubbard U, then addressing W ^> U 
as weak and W%U as strong correlation with W <C U as the atomic limit. Question: How 
are C2,ANx, static/dynamic correlation, W ^ U mutually related? 
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